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Knowledge of the atomic geometry of a surface is a prerequisite for any detailed understanding 
of the surface's electronic structure and chemical properties. Previous studies have convincingly 
demonstrated that density functional theory (DFT) yields accurate surface atomic geometries and 
that reliable predictions concerning stable and metastable phases can be made on the basis of the 
calculated energetics. In the present work we use DFT to investigate the atomic structure of four 
ordered coadsorbate phases of carbon monoxide and oxygen on Ru (0001). All of the structures have 
a (2 x 2) periodicity with differing concentrations of CO molecules and O atoms. For two of these 
phases dynamical low-energy electron diffraction (LEED) intensity analyses have been performed 
and the agreement between our DFT- and the LEED-determined structures is found to be very 
good. We predict the atomic geometry of the third phase for which no structural determination 
based on experiments has been made to date. We also predict the stability of a new ordered mixed 
phase. 



I. INTRODUCTION 

The coadsorption system-|CO,0/Ru (0001) represents 
a well-studied model systcmlllfl for the investigation of the 
interaction and behavior of coadsorbed species that can 
give insight into physicaLprocesses that are relevant to 
heterogeneous catalysisu u. Quite clearly, in order to de- 
scribe and understand chemical reactions at surfaces and 
the surface's electronic structure, it is at first necessary 
to know the atomic geometry. It is only recently that the 
detailed atomic structure of some of the phases of CO,0 
on Ru (0001) have been determined, and new ones discov- 
ered. Depending on the experimental preparation, coad- 
sorption of CO and O on Ru (0001) have been gsported to 
form the following phases: (2 x 2)-(1Q+1C0)BQ, (2 x 2)- 
(2O+1CO)0, and (2 x 2)-(10+2CO)l. In the rest of the 
paper we omit indicating the periodicity since for all of 
the coadsorption structures considered, it is the same. 
For the first two structures, low-energy electron diffrac- 
tion (LEED) intensity analyses have been performed: In 
the first phase, the O atoms occupy hep sites and the CO 
molecule adsorbs in the on-top site; interestingly, a tilt 
of the CO molecular axis of 12.6° was identified. In the 
second phase, a restructuring induced by CO adsorption 
of the O atoms of the (2 x l)u phase occurs: half of the O 
atoms, initially occupying the hep sites, switch to fee sites 
and CO adsorbs again in the on-top site; the other half 
of the O atoms remain in hep sites. For the (10+2CO) 
structure, to our knowledge, there has been no LEED 
intensity analysis, but infrared absorption spectroscopy 
(IRAS) ami X-ray photoelectron spectroscopy (XPS) ex- 
perimentsu indicate that the O atoms occupy hep sites 
and the CO molecules occupy on-top and fee sites. For 
this structure, it was found that in order to achieve the 
higher coverage of CO, notably higher exposures were 
necessary, i.e. ~ 10 5 Langmuirs. In correspondence, the 



sticking coefficient, which is initially close to one, drops 
by a factor of 60 at CO coverages greater than w0.25. 
That is, initially there is a steep CO uptake curve, there- 
after there is a kink and it rises only slowly with exposure. 

Previous studies of the structure and stability of ad- 
sorbate phases using density-functional theory (DFT) 
have demonstrated that atomic geometries and stable 
and metastable structures can be accurately predicted 
(see, e.g£3 Ji_J ) . In the present work we perform DFT 
calculations!! for each of the identified phases in or- 
der to investigate, at first, the atomic structure. This 
work represents one of the first DFT studies of such 
coadsorbate systems where quantitative comparison with 
structures determined by LEED are made. In a subse- 
quent publication we will report on the energetics and 
further details of the electronic structure. We compare 
our calculated atomic geometries for the (10+ ICO) and 
(20+1CO) phases with the results of the LEED inten- 
sity analyses, and predict that of the to date undeter- 
mined (10+2CO) phase. In addition, we investigated a 
hypothetical (30+ ICO) structure, which we predict to 
represent a new high density stable phase. These mixed 
0,CO/Ru (0001) surface structures are depicted in the 
top panel of Fig. 1. 



II. CALCULATION METHOD 

We use the generalized gradient approximation (GGA) 
for the exchange-correlation functionalO-3 and the super- 
cell approach to model the surface structures, which con- 
sists of four Ru layers. The O atoms_and CO molecules 
are adsorbed on one side of the.slablii We use ab initio, 
fully separable pseudopotentialsll3 where the GGA is em- 
ployed for all atomsEO, and relativistic effects taken into 
account for the Ru atoms (using weighted spin-averaged 
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pseudopotentials) . The positions of the C and O atoms, 
and the Ru atoms in the top two layers are relaxed. The 
energy cutoff is 40 Ry and thpe are three special k-points 
in the surface Brillouin zoneli3. 



III. O ON RU (0001) 

In an earlier publication0 we reported results for or- 
dered adlayers of O on Ru(0001); some of these results 
will be used as reference structures for comparisontj. In 
addition, we perform calculations for an artificial "honey- 
comb" (2 x 2)-20 structure in which one O atom in the 
surface unit cell occupies an hep site and the other an 
fee site. The atomic geometry of the honey-comb struc- 
ture is as for the (20+1CO) phase, but without the CO 
molecules. The adsorption of oxygen on Ru (0001) under 
ultra high vacuum-, (UHV) conditions forms two ordered 
phases: (2 x 2)-OE!J and (2 x l)-Oa for coverages 1/4 and 
1/2, respectively. Recently two additional-phases have 
been identified: (1 x l)-OE3 and (2 x 2)-30BB structures 
for coverages 0=1 and 3/4, respectively. Formation of 
the latter two phases require introduction of higher con- 
centrations of oxygen to the surface either by employing 
very high gas pressures of O2 or by using "atomic oxy- 
gen" via NO2 dissociation. In all of the ordered phases 
of O on Ru (0001) that form in nature, O adsorbs in the 
hep site. This is the "natural" site, i.e. the site where 
also Ru would sit; it is stabilized by a lowering of the 
occupied DOS, following the trend that systems like to 
attain a chemically hard electronic structure, i.e. a low 
density at the Fermi level. Thus, essentially the same 
effect which stabilizes-the hep structure for Ru over the 
fee structure (see e.g.f^l). 



three surface Ru atoms. Thus CO adsorption does not 
actually alter the relative position of the Ru atom very 
much. Compared to the vertical distance between the 
center of mass of the first and second Ru layers of the 
(2 x 2)-0 structured, however, CO has induced a slight 
expansion, of 1.2 % as determined by LEED and 1.35 % 
by DFT-GGA. „ 

The LEED intensity analysisu also identified a tilt of 
12.6° of the CO molecular axis in a direction between 
two adsorbed O atoms; it was postulated that the tilt 
was static, rather than dynamic. We investigated such 
a tilting of the CO axis by initially taking the identi- 
fied tilt and allowing atomic relaxation. The final tilt 
angle was found to be 6°. The difference to the LEED 
result lies mainly in the lateral position of the C atom: 
LEED finds a very small displacement (0.05±0.11 A) 
while DFT-GGA obtains a larger value (0.141 A) giving 
rise to a tilt of the C-O-Ru complex rather than just of 
the C-0 axis. The results in any case agree to within the 
experimental error bars. In addition we considered tilting 
of the CO axis towards an adsorbed atom. We took the 
initial tilt to be that of 12.6° as determined by LEED. 
On relaxation of the atomic positions, we found again 
a tilt of «6°. On inspection of the energetics of both 
the tilt-geometries we found that the total energy was 
practically identical to the upright configuration. Thus 
the potential energy surface (PES) is rather flat which 
suggests that the CO molecule will strongly vibrate. 

The corresponding electron density of the valence 
states for this phase is shown in Fig. 1. The oxygen 
atoms appear as the red, almost spherical features. Both 
adsorbate species can be seen to induce a redistribution 
of the electron density of the top-layer Ru atoms, where 
the C-Ru bond is clearly seen. 



IV. (10+1C0)/RU (0001) 

In the coadsorption system, (10+lCCd, CO adsorbs 
in the on-top site and O in an hep siteo. The struc- 
tural parameters determined by LEED and DFT-GGA 
are given in Tab. I (compare Fig. 1). The CO bond 
length determined by both methods is «1.16 A, only 
slightly longer than that of the free molecule which is 
1.15 A (as calculated in the presentpsjork and also as ob- 
tained experimentally, see e.g. Refo) and only-slightly 
shorter than that on the clean surface (wl.17 AEj). The 
LEED-determined C-Ru distance-is 1.93 ±0.06 A (just 
the same as on the clean surfaceEil) and the DFT-GGA 
value is 1.950 A . The position of the Ru atom to which 
CO is adsorbed is displaced slightly inwards by «0.06 A 
(by LEED and DFT-GGA) relative to the tap-most Ru 
atoms in the unit cell. As discussed in Ref.Q, this is in 
apparent contrast to the behavior of CO on the clean sur- 
face where, the Ru atom is "pulled" outwards by 0.07 A . 
We point out, however, that with respect to the (2 x 2)-0 
structure, the same Ru atom to which CO would adsorb 
is already displaced inwards by CUj)5 A as determined by 
LEEDS (0.043 A by DFT-GGAOJ), relative to the other 



V. (20+lCO)/RU (0001) 

The coadsorption system (20+1CO) is formed by CO 
adsorption onto the (2 x l)-0 phase. CO induces a re- 
structuring of the O atoms such that half of them move 
into the fee sites, while the other half remain in hep sites. 
By doing so, CO maintains its favored on-top adsorption 
site. The structural parameters of this phase as deter- 
mined by LEEDQ and DFT-GGA are given in Tab. II 
(compare Fig. 1). It can be seen that there is good 
general agreement. Both LEED and DFT-GGA find a 
slightly shorter CO bond length (i.e., of 1.15 A and 
1.155 A, respectively) than for either CO on the clean 
surface («1.17 A) or in the (10+1CO) phase (ral.16 A). 
There is, however, a slight difference between the LEED 
and DFT-GGA results: The LEED-determined C-Ru 
bond length is longer in this phase than for CO on the 
clean surface or in the (10+1CO) phase (1.98 A com- 
pared to 1.93 A) indicating a weaker C-metal bond- 
ing. The DFT-GGA calculations find a similar, but 
somewhat shorter C-Ru bond length (1.922 A) than in 
the (10+1CO) phase (1.950 A), indicating a somewhat 
stronger bonding of CO to the metal in this phase. In a 
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forthcoming publication we will investigate the energetics 
of all the ordered mixed coadsorption phases. 

The trends in the rumpling and lateral shifts of the top 
layer Ru atoms are agree well between LEED and DFT- 
GGA. Similarly to the (10+1CO) phase discussed above, 
the Ru atom to which CO is bonded is displaced in- 
wards relative to the other top-layer Ru atoms, by 0.09 A 
(LEED) and by 0.04 A (DFT-GGA). For the underlying 
O structure of the (20+ ICO) phase, which corresponds 
to the artificial "honey-comb" geometry described above, 
we point out that a similar, but larger difference already 
exists between the vertical positions of the top-layer Ru 
atoms; namely, the same Ru atom to which CO would 
bond is 0.154 A further in towards the surface than the 
other Ru atoms in the surface unit cell. Therefore, rela- 
tive to the position of the Ru atoms of the honey-comb 
structure, the Ru atom to which CO bonds is in fact dis- 
placed outwards by 0.114 A due to CO adsorption. 

While the LEED study did not report any relaxation of 
the second Ru layer, the DFT-GGA calculations identify 
small relaxations and a rumpling: The Ru atom below 
the hep O atom is displaced 0.038 A further in towards 
the bulk than the other three Ru atoms. These latter 
three Ru atoms exhibit lateral displacements of 0.013 A 
radially in towards each other. The top Ru-Ru verti- 
cal displacement, using the center of mass of each layer, 
is 0.6 % larger than that of the honey-comb geometry, 
which is due to CO adsorption. 

From the valence electron density (Fig. 1) it can be 
seen that the bonding of the hep O atom to the metal 
appears to be stronger than the fee O atom; this is what 
we may expect on the basis of the magnitude of the en- 
ergy difference of between the hep and fee sites in the 
stable (2 x 1) structure for the same coverageEd. 

VI. (lO+2CO)/RU(0001) 

For the (10+2CO) phase no LEED analysis exists 
so far. The calculated atomic geometry is depicted in 
Fig. 1. The CO bond length for CO in the on-top site is 
1.157 A and for CO in the fee site the bond length is a 
longer 1.187 A. For the top site, the C-Ru bond length 
is 1.968 A, similar to that of the (10+1CO) phase (and 
slightly longer than that in the (20+1CO) phase, which 
was 1.922 A as determined by DFT-GGA). For the fee 
site the C-Ru bond length is 2.238 A , notably longer 
than for CO in the on-top site; this, however, is expected 
due to the fact that it is three-fold coordinated and not 
one-fold. The calculations also identify small lateral dis- 
placements of 0.018 A of the three top-layer Ru atoms 
in the direction towards the fee site. Also, a rumpling of 
the top Ru layer of 0.129 A occurs as indicated in Fig. 1, 
where again the Ru atom to which CO bonds in the on- 
top site is further in towards the bulk. In this case, with 
respect to the (2 x 2)-0 structure, the Ru atom is dis- 
placed inwards by 0.086 A (i.e., 0.129-0.043 A ). In the 
second Ru layer, the Ru atoms below the hep O atoms 
are 0.021 A in further towards the bulk than the other 



three Ru atoms, and again these three Ru atoms dis- 
play a lateral radial displacement towards each other of 
0.023 A . We hope this predicted atomic geometry can be 
confirmed by a LEED intensity analysis. 

From the cross-section shown of the valence electron 
density, the bonding of the CO molecule in the fee site 
can be seen to be weaker than for that of the on-top site. 
In the fee site, however, CO forms three bonds with the 
metal surface so it is expected that they be longer and 
weaker. 



VII. (30+lCO)/RU (0001) 

We have also studied the atomic geometry of a 
(30+1CO) structure, although, to date, such a phase 
has not been reported to exist experimentally. Our re- 
sults predict that the structure is energetically stable, 
although the adsorption energy is low: Starting from a 
(2 x 2) array of O vacancies in the full oxygen mono- 
layer the CO adsorption energy is 0.85 eV. As mentioned 
above, recently it has been shown that such a (2 x 2)- 
3O/Ru(0001) "vacancy" or "hole" structure ip-£act rep- 
resents a new stable phase of O on Ru (0001)Ei}El The 
stability of this phase was also indicated in our studyj-of 
the thermodynamics of the O/Ru adsorption systcmc3. 
In order to investigate the energetics of CO adsorption 
into the vacant hep site we calculated the total energy of 
CO at various distances above this site. We found that 
there is an energy barrier of ~0.3 eV to adsorption. This 
implies that rather high CO pressures would be required 
in order to realize this structure experimentally. Provid- 
ing these activation barriers can be overcome, we propose 
that this structure represents a new ordered (2 x 2) mixed 
phase of CO and O on Ru (0001). The atomic geometry is 
depicted in Fig. 1. Here the CO bond length is 1.181 A, 
somewhat longer than for CO in the on-top site, but sim- 
ilar to that of CO in the fee site for the previous structure 
discussed, which was 1.187 A . In the (30+1CO) system, 
there are only negligible lateral shifts of the O atoms and 
also only negligible rumpling of the top Ru layer (i.e., all 
<0.005 A). In the second Ru layer, however, the rum- 
pling is 0.043 A . In this case the Ru atom directly under 
CO is displaced outwards (or upwards) relative to the 
other three Ru atoms. There are again only negligible 
lateral displacements in the second Ru layer. Due to CO 
adsorption, using the center of gravity of the layers, the 
vertical Ru-Ru distance is expanded by «2 % relative 
to the (2 x 2)-30 O-vacancy structure. Thus again, CO 
adsorption has induced an overall expansion of the first 
Ru-Ru interlayer spacing. 

The valence electron density in Fig. 1 shows that there 
is the most disturbance of the Ru states for this structure. 
With respect to CO adsorption, we see four clear maxima 
of the electron density on the Ru atom below, while for 
the Ru atom that only forms bonds with O atoms, two 
strong maxima appear which are polarized away from the 
O-Ru bond. This redistribution of Ru states is seen to 
give rise to a build-up of electron density on top of the 
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Ru atoms in the top layer. 



VIII. CONCLUSION 



We have investigated mixed O, CO adlayers on 
Ru(0001). We find very good overall agreement of 
the atomic geometry of the (10+1CO) and (20+1CO) 
phases with those determined by LEED intensity anal- 
yses. We predict the to date undetermined structural 
geometry of the (10+2CO) phase and propose the exis- 
tence of a new high density phase, namely, (30+1CO). 
With respect to the latter structure, adsorption of CO 
onto the (2 x 2)-30 surface in order to achieve this phase 
is found to be activated. This, together with the notably 
reduced area for CO adsorption on this surface, suggests 
that high CO pressures would be required to create this 
phase experimentally. 



4 



(nO+mCO) Adsorbate Phases on Ru (0001) with (2x2) Periodicity 



(lO+ICO) 



(20+ ICO) 





(10+2CO) 





d(O-RuII) 

d(O-RuII) 



z(Ru A -Ru) 




FIG. 1. Perspective and side views of the various phases of O and CO on Ru (0001). Large and small (green and red) circles 
represent Ru, O, and C atoms, respectively. The lower panel shows the electron density of the valence states. The contour 
lines are given in e bohr -3 . 





d(C-O) 


d(C-Ru) 


d(0-Ru A ) 


d(0-Ru B ) 


z(Rua-Ru) 


A||(Ru A ) 


A,, (Ru B ) 




(A) 


(A) 


(A) 


(A) 


' (A) 


(A) 


(A) ' 


DFT-GGA 


1.161 


1.950 


2.109 


2.166 


0.057 


0.050 


0.068 


LEED 


1.16±0.06 


1.93±0.06 


2.06±0.08 


2.09±0.14 


0.06±0.05 


0.16±0.10 


0.10±0.10 



TABLE I. Structural parameters obtained by DFT-GGA and by LEEDu for the (lO+ICO) phase (see Fig. 1). 





d(C-O) 


d(C-Ru) 


z(0(fcc)-0(hcp)) 


d(0(hcp)-RuII) 


d(0(fcc)-RuII) 


z(RuII-RuI) 


Ay (Run) 




(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


DFT-GGA 


1.155 


1.922 


0.157 


2.080 


2.122 


0.040 


0.039 


LEED 


1.15±0.04 


1.98±0.08 


0.19±0.09 


2.05±0.16 


2.05±0.16 


0.09±0.02 


0.08±0.08 



TABLE II. Structural parameters obtained by DFT-GGA and by LEED0 for the (20+1CO) phase (see Fig. 1). 
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